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Abstract 

We present "state-of-the-art" theoretical expressions for the triple differential B — » 
X u l~v decay rate and for the B — > X s ^f photon spectrum, which incorporate all known 
contributions and smoothly interpolate between the "shape-function region" of large 
hadronic energy and small invariant mass, and the "OPE region" in which all hadronic 
kinematical variables scale with Mb- The differential rates are given in a form which 
has no explicit reference to the mass of the b quark, avoiding the associated uncertain- 
ties. Dependence on m& enters indirectly through the properties of the leading shape 
function, which can be determined by fitting the B — > X s ^ photon spectrum. This 
eliminates the dominant theoretical uncertainties from predictions for B — > X u l~v de- 
cay distributions, allowing for a precise determination of \V u b\. In the shape- function 
region, short-distance and long-distance contributions are factorized at next-to-leading 
order in renormalization-group improved perturbation theory. Higher-order power cor- 
rections include effects from subleading shape functions where they are known. When 
integrated over sufficiently large portions in phase space, our results reduce to standard 
OPE expressions up to yet unknown 0{ot 2 s ) terms. Predictions are presented for partial 
B — > X u l~v decay rates with various experimental cuts. An elaborate error analysis is 
performed that contains all significant theoretical uncertainties, including weak annihi- 
lation effects. We suggest that the latter can be eliminated by imposing a cut on high 
leptonic invariant mass. 



1 Introduction 



A major effort of the 5-physics community is underway to map out the apex of the unitarity 
triangle, which provides a graphical representation of the effect of CP violation in the quark 
flavor sector of the Standard Model. One of the biggest successes of this endeavor was the 
precise determination of the angle /3, which has been measured with high accuracy from the 
time-dependent CP asymmetry in the B — > J/ipKs decay channel [1, 2]. The length of the side 
opposite the angle j3 is proportional to |14&|- A high-precision determination of this quantity 
would enable us to test the validity of the Standard Model and search for possible deviations 
from its predictions. 

Good theoretical knowledge of strong-interaction effects in weak decays of B mesons is 
crucial for a reliable exploration of the flavor sector of the Standard Model. In particular, 
the determination of the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements |V^,| and |T4&| 
relies on an accurate description of bound-state effects in semileptonic decays. At present, the 
most precise calculations are available for inclusive semileptonic decays B — > Xl~u. 

The theoretical tools for the calculation of inclusive B decays are QCD factorization on 
the one hand [3, 4, 5, 6, 7, 8, 9, 10, 11], and local operator product expansions (OPE) on the 
other [12, 13]. Both approaches perform a systematic separation of long-distance hadronic 
quantities from short-distance perturbative ones, while organizing the calculation in inverse 
powers of the heavy 6-quark mass m&. The OPE is an appropriate tool for the calculation 
of total inclusive rates (for example in B — > X c l~u decay) or for partial rates integrated 
over sufficiently large regions in phase space, where all components of the final-state hadronic 
momentum P x are large compared to Aq CD . QCD factorization, on the other hand, is better 
suited for the calculation of partial rates and spectra near kinematical boundaries, where 
typically some components of P x are large, while the invariant hadronic mass M x = a/ P x 
is small. For example, any B — > X u l~v event can be described with three independent 
kinematical variables, a useful choice of which is [9, 14] 

P l = M B -2E l , P_ = E X + \P X \, P + = E X -\P X \. (1) 

Here P± are the light-cone components of the hadronic final-state momentum along the jet 
direction, E x is the charged-lepton energy, E x is the jet energy, and P x is the jet momentum, 
all measured in the 5-meson rest frame. The phase space for these variables is 

M 2 

-f- < P+ < Pi < P- < M B , (2) 

with M w being the mass of the lightest possible hadronic final state. The product P + P- = M\ 
is the hadronic invariant mass squared. In order to avoid large backgrounds from b — > c 
transitions, all measurements of \V u b\ are in one way or another restricted to the region of 
phase space where P+P- < M^. If the quantity P_ takes values near its maximum at M B , 
then P + is restricted to a region of order Mjj/M b , which is numerically comparable to A QCD . 
This means that there are three parametrically different energy scales in the problem: the mass 
M B of the initial state, the mass of the final hadronic state ~ MbAqcd 5 and the low scale 
Aqcd at which perturbation theory breaks down and hadronic physics must be parameterized 
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in terms of non-perturbative matrix elements. QCD factorization disentangles the effects from 
these scales, so that the perturbative contributions can be expanded in powers of a s (nh) with 
jih ~ m h (giving rise to "hard functions") and a s (//j) with /ij ~ a/"V^qcd (giving rise to "jet 
functions"). 

It is important to note that the heavy-quark expansions valid in these two kinematical 
regions are not identical, because the power counting rules differ in the two regimes. Also the 
nature of the non-perturbative inputs is different. In the OPE region, non-perturbative physics 
is encoded in a few hadronic parameters, and the heavy-quark expansion is the usual Wilsonian 
expansion in local operators. In the endpoint (or shape-function) region, the presence of 
multiple scales complicates the power counting, and the interplay between soft and collinear 
modes gives rise to large non-localities. As a result, non-perturbative physics is described by 
hadronic structure functions called "shape functions", and the heavy-quark expansion is an 
expansion in non-local string operators defined on the light-cone. The connections between 
the two regimes is that moments of the shape functions can be expressed in terms of local 
operators. 

The goal of the present work is to develop a formalism that smoothly interpolates between 
the two kinematical regimes (see [15] for a related discussion, which is however restricted 
to the tree approximation). This is essential for building an event generator for inclusive 
B — > X u l~v and B — > X s ^ decays, which can be used to study partial and differential decay 
rates in different kinematical domains. In the shape-function region, our approach relies on 
exact QCD factorization theorems, which exist in every order of power counting. They allow 
us to systematically disentangle short- and long-distance physics and, in the process, resum 
parametrically large logarithms order by order in perturbation theory. This factorization can 
be done with high accuracy for the terms of leading power in l/m&, and with somewhat less 
sophistication for the first-order power corrections. For the second-order power corrections, we 
only include contributions that do not vanish when integrated over all phase space. This is a 
safe approximation; the effects of the remaining l/m% terms can to a large extent be absorbed 
by a redefinition of the subleading shape functions arising at order l/m^. 

Our formalism is "optimized" for the shape-function region in the sense that sophisticated 
theoretical technology is applied in this regime. However, when our expressions for the differ- 
ential decay rates are integrated over sufficiently wide domains, they automatically reduce to 
the simpler results that can be derived using the OPE approach, up to yet unknown terms of 
0(a 2 s ). The moment relations for the shape functions are crucial in this context. Note that 
local 1/ml corrections in the OPE receive contributions from terms of leading power (l/m°), 
subleading power (l/m^), and sub-subleading power (1/mf) in the shape-function region, so 
the transition is highly non-trivial. In implementing the program outlined here, we include all 
presently known information on the triple differential B — > X u l~v decay rate and on the differ- 
ential B — > X s 7 decay rate in a single, unified framework. We neglect, for simplicity, hadronic 
power corrections of order l/mf and higher, which are known to have a negligible effect on 
the observables considered here. The only possible exception is contributions from "weak 
annihilation", which are estimated as part of our error analysis. We also ignore the existing 
results on O(f3oa 2 ) radiative corrections for some single-differential distributions, because the 
corresponding corrections are not known for the double or triple differential B — > X u l~v decay 
spectra. While these O(P al) terms are sometimes found to be large when naive perturbation 
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theory in a s {m b ) is used, their effects are expected to be small in our scheme, which is based 
on a complete scale separation using QCD factorization. We see no reason why the floa^ terms 
should be enhanced compared to other, unknown corrections of 0(af). 

A technical complication in realizing the approach described here has to do with the 
treatment of phase-space factors. The heavy-quark expansion of the hadronic tensor for B — > 
X u l~u decay gives rise to expressions that are singular at certain points in phase space. One 
way to avoid these singularities is to also expand phase-space factors order by order in Xjm^ 
(see, e.g., the treatment in [16]). However, since this expansion depends on the kinematical cuts 
of any given analysis, it cannot be implemented in a straightforward way in an event generator. 
An alternative is to reorganize the heavy-quark expansion in such a way that the expansion 
parameter is related to hadronic (as opposed to partonic) kinematical variables, in which case 
kinematical singularities are always canceled by exact phase-space factors. Following this 
strategy, we obtain expressions for decay distributions and partial decay rates which are free 
of explicit reference to partonic quantities such as the 6-quark mass. A dependence on to& 
enters only implicitly via the first moment of the leading-order shape function S{CS). The 
philosophy of our approach is that this function 1 is extracted experimentally from a fit to the 
B — > X s 7 photon spectrum, which has been measured with good precision in the region where 
P + = Mb — ~ Aqcd- This is analogous to the extraction of parton distribution functions 
from deep inelastic scattering. The photon spectrum is experimentally accessible to energies 
as low as 1.8 GeV, which corresponds to a sampling of the shape function for values of Co up 
to about 1.7 GeV. Once the shape function has been extracted over this range, we can use it 
to obtain predictions for arbitrary partial B — > X u l~v decay rates with cuts. In doing so, 
the residual hadronic uncertainties in the extraction of \V u b\ only enter at the level of power 
corrections. 

We emphasize that the program outlined above is equivalent to an approach put forward 
in [4] and later refined in [17, 18, 19], in which \V u b\ is extracted with the help of shape- 
function independent relations between weighted integrals over differential decay distributions 
in B — > Xg'-f and B — > X u l~v. The experimental error in the results for these weighted 
integrals corresponds, in our approach, to the error in the prediction of B — > X u l~v partial 
rates resulting from the experimental uncertainty in the extraction of the shape function from 
the B — > X s 7 photon spectrum. While the shape-function independent relations are very 
elegant, it is more convenient for the construction of a generator to have a formulation where 
the shape function is used as an input. In this way, it is possible to impose arbitrary cuts on 
kinematical variables without having to recompute the weight functions in each case. 

The paper is structured as follows: In Section 2 we collect the relevant formulae for the 
calculation of the B — > X s ^y photon spectrum. These expressions can be used to extract the 
leading non-perturbative structure function from experiment. An analogous presentation for 
the triple differential decay rate in B — > X u l~v decays is presented in Section 3. In order to 
perform a numerical analysis one needs to rely on parameterizations of the shape functions. A 
collection of several useful functional forms is given in Section 4. In Section 5 we present a full 

1 More precisely, we define a new shape function S(uj) by the combination of leading and subleading shape 
functions contributing to B — > X s j decay, and we will use the same function to make predictions for B — > 
X u l~v decay distributions. 
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error analysis of partial B — > X u l~u decay rates for a variety of experimental cuts. We also 
explore the sensitivity of the results to the 5-quark mass and to the functional forms adopted 
for the shape functions. Section 6 contains our conclusions. 

2 Inclusive radiative decays 

The decay process B — > X s 7, while more complex in its short-distance physics, is considerably 
simpler in its kinematics than the semileptonic process B — > X u l~u. Since the radiated 
photon is on-shell, the hadronic variables P± that describe the momentum of the X s system 
are trivially related to the photon energy E 1 by P + = Mb — 2E^ and P_ = Mb- In the crudest 
approximation, namely at tree level and leading power, the photon-energy spectrum is directly 
proportional to the leading shape function, dT s /dE^ oc S(P + ). In this section we collect all 
relevant formulae needed to compute the B — > X s ^ photon spectrum or, equivalently, the 
invariant hadronic mass distribution. It is implicitly assumed that these spectra are sufficiently 
"smeared" (e.g., by experimental resolution) to wash out any sharp hadronic structures. In 
cases where the resolution is such that the K* resonance peak is observed, it can be accounted 
for by combining the formulae in this section with the prescription for subtracting the K* 
peak proposed in [20] . 

The differential B — > X s 7 decay rate can be written as 

= ^|? e 7 3 iw:i 2 ^m [ c ^(t i h)] 2 u(n h ,tii)^(p + ), (3) 

where the structure function jF 7 depends on the photon energy via P + = Mb — 2E 1 . The pref- 
actor contains the electromagnetic fine-structure constant a normalized at q 2 = 0, two powers 
of the running 6-quark mass (defined in the MS scheme) originating from the electromagnetic 
dipole operator Q 7l in the effective weak Hamiltonian, and the square of the corresponding 
Wilson coefficient C^, which is needed at next-to-leading order in renormalizat ion-group im- 
proved perturbation theory [21]. Renormalization-group running from the hard scale fih ~ n^b 
to the intermediate scale /ij ~ a/"V^cd gives rise to the evolution factor U (fih, Hi), whose 
explicit form is discussed in Appendix A. We keep U and (C^) 2 outside of the structure 
function HF 7 ; it is understood that when combining the various terms in (3) all perturbative 
quantities should be expanded for consistency to the required order in a s . 

2.1 Leading-power factorization formula 

At leading order in 1/mj, the structure function jF 7 factorizes as [11] 

rP+ 

3f{P + ) = \H s (^ h )\ 2 / dujm h J{m h {P + -uj), i i i )S{Cj^ i ). (4) 

Jo 

At this order a single non-perturbative parton distribution function arises, called the leading 
shape function [4] and denoted by S(cu, /ij). Our notation is adopted from [9, 16]: hatted shape 
functions have support for Co > 0. The function S is defined in terms of a non-local matrix 
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element in heavy-quark effective theory (HQET). Renormalization-group running between 
the intermediate scale and a low hadronic scale is avoided when using the shape functions 
renormalized at the intermediate scale /ij. Evolution effects below this scale are universal (i.e., 
process independent) and so can be absorbed into the renormalized shape function. Short- 
distance contributions from scales above \ih ~ m& are included in the hard function H s , which 
in practice is obtained by matching the effective weak Hamiltonian onto a current operator 
in soft-collinear effective theory (SCET). At next-to-leading order in perturbation theory, the 
result reads 



H s (n h ) = 1 + 



C F a s (n h 



-2 In 2 



+ 



+ 



4tt \ fi h 
c t g {^h) C F a s (fi h ) 



47T 



- + 71n 6-- 



2ni 



8, m b 11 2tH 

— In 1 h 

3 jih 3 9 3 



r "" y -[9(0)-9W] 



104 , , , „ 



v tb v t * s 



+ e 



pcng 



, (5) 



where the variable z = (m c /mb) 2 denotes the ratio of quark masses relevant to charm-loop 
penguin diagrams, and the "penguin function" g(z) can be approximated by the first few 
terms of its Taylor expansion, 



833 20m 8tt 2 
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,3/2 
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162 27 

+ _£ [48 - 5tt 2 - 36Ca + (30tt - 2n 3 )i + (36 - 9tt 2 + 6m) In z + (3 + 6m) In 2 z + In 3 z] 
9 

2z 2 

+ — [18 + 2tt 2 - 2n 3 i + (12 - 6tt 2 ) In z + 67ri In 2 * + In 3 z] 
9 



+ [-9 -14tt 2 + 112^ + (182 -48^) In 2 -126 In 2 2] + ... . 

— I 



(6) 



The Wilson coefficients C\ and C|J in (5) multiply the current-current operators Qi ,c and 
the chromo-magnetic dipole operator Q 8g in the effective weak Hamiltonian. The quantities 



■1.5% and e 



peng 



-0.6% account for small electroweak corrections and the effects of 
penguin contractions of operators other than Q"' c , respectively. The differential decay rate (3) 
is formally independent of the matching scales fih and /ij. The fih dependence of the evolution 
factor U(nh,t^i) cancels the scale dependence of the product m 2 (/i^) [Cj^(fih)} 2 \H s (nh) \ 2 , while 
its fii dependence compensates the scale dependence of the convolution integral J(//j) <S>S(ni). 

Finally let us discuss the jet function J, which appears as the hard-scattering kernel in the 
convolution integral in (4). It can be written in terms of distributions that act on the shape 
function S. At one-loop order, the jet function is given by [8, 9] 



J(p 2 ,/i) = % 2 ) 



C F a s {fj.) 2 
1 + ~ i^* 7 " 7 ^ 



+ 



C F a s (fx) 

47T 



1 / p 2 



M^ 2 



(7) 
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where the star distributions have the following effect on a function / when integrated over a 
domain Q 2 [22]: 



dp 2 



<o 



rQ 2 

/ dp 2 

J<0 



1 f 

P 2 jjL z 



\M 2 



y 2 



+ /(0) In 



f(p z ) = / dp 

'0 



Q2 ,2 f(P 2 ) - /(0) , /(0) , , Q 



p 1 



ln^ + 



ln^ 



(8) 



2.2 Kinematical power corrections 

There exists a class of power corrections to (4) that do not involve new hadronic quantities. 
Instead, the power suppression results from the restriction of certain variables (P + in the 
present case) to a region where they are kinematically suppressed (here P + <C M B ). The 
corresponding terms are known in fixed-order perturbation theory, without scale separation 
and renormalization-group resummation [23, 24] (see also [20]). To perform a complete RG 
analysis of even the first-order terms in 1/ m& is beyond the scope of the present work. Since, 
as we will see later, power corrections only account for small corrections to the decay rates, 
an approximate treatment will suffice. To motivate it, we note the following two facts [11]: 
First, while the anomalous dimensions of the relevant subleading SCET and HQET operators 
are only known for a few cases [25], the leading Sudakov double logarithms are the same as 
for the terms of leading power, because they have a geometric origin in terms of Wilson lines 
[26]. The leading Sudakov double logarithms are therefore the same as those resummed into 
the function U in (3). Secondly, the kinematical power corrections in B — > X s ^ decay are 
associated with gluon emission into the hadronic final state X s . Because of the kinematical 
restriction to low-mass final states, i.e. M\ ~ M b Aq CD , we associate a coupling a 8 {fi) with 
these terms, where typically fi ~ //j. Strictly speaking, however, the scale ambiguity associated 
with the choice of p, could only be resolved by computing the relevant anomalous dimensions. 

Within this approximation, the kinematical power corrections to the structure function 5F 7 
can be extracted from [11, 20]. We find it convenient to express the result in terms of the 
variable 

which in the shape-function region scales like Aqcd/^6- We obtain 



1 C F a s {n) 



M R -P+ 



An 



E 



A2 Ci{n h ) 



i,j = 1,7, 
i < j 



duo S(cu, hij{x) 



9m 2 c CfJn h ) 



S{P+,i*). 



(10) 
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The coefficient functions hij(x) are 
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where as before z = (m c /mb) 2 , and 



-2arctanVt/(4-t) 



z(l + x)J 2 
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(11) 



G(t) = 



2 In 



(12) 



In the shape-function region the expressions for hij(x) could, if desired, be expanded in a 
power series in x = 0(AQ CD /m b ), and this would generate a series of power-suppressed terms 
g-km(n)^p^ with n > 1, where the superscript "n" indicates the order in the l/m& expansion. 
Note that this expansion would contain single logarithms In a; ~ I^Aqcd/^)- These are 
precisely the logarithms that would be resummed in a more proper treatment using effective 
field-theory methods. 

Outside the shape-function region the variable x can take on arbitrarily large positive 
values, and 3^ in (P + ) is no longer power suppressed. Note that for P + — > Mb (corresponding 
to x — > oo and E 1 — > 0) most functions hij(x) grow like x 2 or weaker, so that the spectrum 
tends to a constant. The only (well known) exception is hss(x), which grows like a; 3 , giving rise 
to a 1/Ery soft-photon singularity [24]. The main effect of the kinematical power corrections 
(10) to the photon spectrum is to add a radiative tail extending into the region of small photon 
energies. These corrections therefore become the more significant the larger the integration 
domain over .EL is. 



2.3 Subleading shape- function contributions 

At order l/nib in power counting, different combinations of subleading shape functions en- 
ter the B — > X s 7 and B — > X u l~v decay distributions [27, 28, 29, 30]. They provide the 
dominant hadronic power corrections, which must be combined with the kinematical power 
corrections discussed in the previous section. We include their effects using the results of 
recent calculations in [16, 31, 32]. Little is known about the subleading shape functions apart 
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from expressions for their first few moments. In particular, the norms of these functions van- 
ish at tree level, while their first moments are determined by the HQET parameters Ai and 
A2, which are defined via the forward S-meson matrix elements of the kinetic-energy and the 
chromo- magnetic operators, respectively [33]. 

For the case of B — > X s j decay, subleading shape-function contributions are currently 
only known for the matrix elements of the dipole operator Q 7 ^, and the corresponding hard 
and jet functions have been computed at tree level. Adopting the notations of [16], the 
relevant subleading shape functions are £(u>), and v{uj). An additional function, called 

So, has been absorbed by a redefinition of the leading shape function, and it is included in 
our definition of S(uj). Roughly speaking, u(uj) is the "light-cone generalization" of the local 
HQET kinetic-energy operator. The functions v(Co) and i(uo) are both generalizations of the 
local chromo- magnetic HQET operator, but t(uj) contains also a light-cone chro mo-electric 
operator, which has no equivalent in the local OPE expansion. (Such a contribution arises 
since there are two external 4-vectors in the SCET expansion, n and v, while there is only v 
in the HQET expansion.) The contribution of subleading shape functions to the B — > X s ^y 
photon spectrum is 

3*-*<i) (P+) = __L_ [_(A - P + ) 5(P + ) -t(P + ) +u(P + ) -v(P + )] . (13) 

Compared to [16], we have replaced 1/mb with 1/(M#— P + ) in the prefactor, which is legitimate 
at this order. (The form of the shape functions restricts P + to be of order Aq CD .) The 
appearance of the HQET parameter A = (M B — m^m^oo is peculiar to the subleading shape- 
function contributions. This quantity is defined via the first moment of the leading-order 
shape function [9]. 

The formula given above can be modified to suit the purpose of extracting the shape 
function from the photon spectrum better. To this end, we absorb a linear combination of 
the subleading shape functions into a redefinition of the leading shape function, in such a way 
that the moment relations for this function remain unchanged to the order we are working. 
This is accomplished by defining 

msm+ mz*m=m±m=m. (14) 

m b 

When using S instead of S in the leading-power formula (4), the subleading shape-function 
contribution becomes 

^) {P+) = J_^_lAs{P + ). (15) 

The hatted shape functions used in the present work are related to the original definitions 
in [16] by 

oiuj) = S(A — lj) H , 

i(u>) = t(A - u) , u(u) = u(A - Cj) , v(u>) = v(A - u) , (16) 
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where the unhatted functions have support on the interval between — oo and A. It is convenient 
to rewrite A — uj = uj + Auj, where 

A -4- 3 A 

Au = A — (Mb — m b ) = 1 J~ 2 + ... (17) 

2m b 

accounts for the mismatch between the HQET parameter A and the difference (Mb — mb) 
due to power-suppressed terms in the l/ra& expansion [34]. It follows that the variable uj = 
(Mb — rrib) — uj runs from — oo to (Mb — Tib)- The moment relations for the leading and 
subleading shape functions derived in [4] and [16, 27] can be summarized as 

ni a \ so(uj + Auj) Ai .... . Ai + 3A2 s 

S(uj) = S(uj + Auj) + '- = 5(uj) - -± 5"(uj) + S'(u) + . . . , 

mb 6 Zrrib 

i(uj) = t(uj + Auj) = A 2 5'(uj) + ... , 

2Ai 

u(uj) = u(uj + Auj) = — S'(uj) + . . . , 

o 

v(uj) = v(uj + Auj) = -A 2 5'(uj) + . . . . (18) 

The function <S has the same moment expansion as S. The hadronic parameter A 2 determines 
the leading contribution to the hyperfine splitting between the masses of B and B* mesons 
through m%, - m% = 4A 2 + 0(l/m b ) [33], from which it follows that A 2 « 0.12 GeV 2 . The 
value of the parameter Ai is more uncertain. In much the same way as the 6-quark pole mass, 
it is affected by infrared renormalon ambiguities [35, 36]. It is therefore better to eliminate Ai 
in favor of some observable, for which we will choose the width of the leading shape function. 



2.4 Residual hadronic power corrections 

At order 1 / ml a new set of sub-subleading shape functions enter, which so far have not been 
classified completely in the literature. Since the functional form of even the subleading shape 
functions is rather uncertain, there is no need to worry too much about the precise form of 
sub-subleading shape functions. Most of their effects can be absorbed into the subleading 
functions. An exception, however, are terms that survive when the sub-subleading shape 
functions are integrated over a wide domain. Whereas the norms of all subleading (~ 1/mb) 
shape functions vanish, the norms of the sub-subleading shape functions (~ l/ m b) are m 
general non-zero and given in terms of the heavy-quark parameters Ai and A 2 . (At tree level, 
the class of functions with non-zero norm has been studied in [15].) Our strategy in the 
present work will be as follows: We start from the well-known expressions for the (tree-level) 
second-order power corrections to the B — > X s 7 photon spectrum [37] (and similarly for the 
triple-differential B — > X u l~v decay distribution [12, 13], see Section 3.4). They are of the 
form \ jm\ times one of the singular distributions S(p 2 ), S'(p 2 ), or S"(p 2 ), where p 2 = (m b v—q) 2 
is the invariant partonic mass squared of the final-state jet. As mentioned earlier, the power 
counting in the shape-function region is different from the one used in OPE calculations, 
and indeed a good portion of the l/m 2 . terms in the OPE is already accounted for by the 
contributions proportional to the leading and subleading shape functions in (4) and (13). We 
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identify the corresponding terms using the moment relations for the shape functions in (18). 
In particular, this reproduces all terms at order 1/mf in the OPE which contain derivatives 
of S(p 2 ). We include the remaining terms of the form (Xi/m 2 ) S(p 2 ) by replacing 



5(p 2 ) = 8{p + p-) = [ du5(p + + u)5(u) 

P--P+ J 



_ , dw8{P + -u)S(u) = J*( P+ l . (19) 
p_- p + J y + 1 y 1 P_ - P + y ' 

Here p± are the light-cone projections of the partonic momentum p^ 1 , which are related to the 
hadronic quantities P± by P± = p± + (Mb — mj,). Similarly, Co = (Mb — rrib) — ou. 
The result of these manipulations is 

^ rW= (M B -P + )^ (P+) - (20) 

Together with (4) and (13) this accounts for all known first- and second-order power corrections 
to the B — > X s 7 photon spectrum, both in the shape-function region and in the OPE region. 
The redefinition (14) of the leading shape function from S to S leaves the form of the second- 
order power corrections unaffected. 

In Section 5 we study the numerical impact of second-order power corrections on various 
B — > X u l~v partial rates and find their effects to be tiny. It is therefore a safe approximation 
to neglect hadronic power corrections of order l/mf or higher. The only possible exception to 
this conclusion relates to the so-called weak annihilation terms in B — > X u l~u decay, which 
will be included in our error analysis. 



3 Inclusive semileptonic decays 

All hadronic physics in B — > X u l~v decays is encoded in the hadronic tensor W^, which is 
defined via the discontinuity of the forward 5-meson matrix element of a correlator of two 
flavor- changing weak currents J M = M7 M (1 — 75) b. Explicitly, 

W^ v = — Im (B(v)\i [ d 4 xe iq x T{j^(0),r(x)}\B(v)) , (21) 
2Mb tt J 

where v is the 5-meson velocity and q the momentum carried by the lepton pair. The hadronic 
tensor can be decomposed into five structure functions Wi, which are the coefficients of the 
five possible Lorentz structures built out of two independent 4-vectors. Typical choices for 
these two vectors are q and v, p and v, etc. Here, as above, p = m^v — q is the momentum of 
the jet of light particles into which the b quark decays. In principle, all choices are equivalent, 
and it is solely a matter of convenience which basis one picks. 

The triple differential decay rate can then be expressed in terms of kinematical prefactors 
and the functions Wi. It is a known fact that the total decay rate is proportional to five 
powers of the 6-quark mass. Further sensitivity to m& is picked up for partial decay rates by 
the kinematical cuts. For example, cutting on the leptonic invariant mass q 2 > q$ introduces 
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roughly five additional powers, and the resulting partial decay rate is proportional to {m^Y 
with a ^ 10 [38, 39]. This is the reason why theoretical predictions were typically made 
for event fractions, so that at least the five powers of mj in the total rate drop out. For 
practical purposes, however, this procedure presents no advantage as the value of the total 
decay rate cannot be measured. Furthermore, the mj dependence of the total rate is clearly 
related to the ra^ dependence of partial rates, and it is important to take this correlation into 
account when combining calculations of event fractions with those of the total decay rate. In 
Section 5, where we present theoretical predictions, we will thus focus directly on predictions 
for partial decay rates, not event fractions. Note that information about enters the triple 
differential decay rate in two ways, via the hadronic structure functions Wi and through their 
kinematical prefactors. Whether or not m b appears explicitly in the prefactors depends on 
the decomposition of W^ u , i.e., on the choice of vectors used to form the five possible Lorentz 
structures. 

A very useful set of 4- vectors turns out to be (v,n), where n is a light-like vector in the 
direction of the jet of light particles. In SCET, n denotes the direction of the collinear particles 
in the jet, which is typically set to be along the z-axis. The normalization is chosen such that 
v ■ n — 1, so that n M = (1, 0, 0, 1) in the rest frame of the B meson. The conjugate direction to 
n is denoted by n M = (1, 0, 0, —1) and marks the direction of the photon in B — > X s 7 decay, 
or the direction of the lepton pair in B — > X u l~~v decay. We then decompose 

= {itv v + n v xP - - ie^ al3 n aV(3 ) W x - g^ W 2 

+ u'V W 3 + {n»v v + n 1 ^) W 4 + n»n v W 5 . (22) 

The structure functions Wi all have mass dimension —1 in this basis. In terms of the Wi 
functions the triple differential decay rate reads 



ub 



- U y (fi h , (M B - P + ) (P_ - P t )(M B - P_ + P l - P + ) ^ 



dP + dP^ dPi 16tt 3 

+ (M B - P-)(P- - P+) ?2 4 (P- - Pi)(Pi - P+) ? 3 ] , (23) 

where we have collected the relevant combinations of Wi into the three functions 

(24) 

and defined a new kinematical variable 

» = 5^%- (25) 

which can take values < y < 1. The leading evolution factor Uy(/j,h, to) has been factored 
out in (23) for convenience, as we have done earlier in (3). The function U y (toi, to) differs from 
the corresponding function in B — > X s ^ decay by a y-dependent factor, 

U v (» h , to) = U(to, to) y- 2a ^ M , (26) 
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where the function ar in the exponent is related to the cusp anomalous dimension and is given 
in Appendix A. 

Eq. (23) for the triple differential rate is exact. Note that there is no reference to the 
6-quark mass at this point. The only dependence on m b is through the structure functions 
$i(P + , y) (via hard matching corrections and via the moment constraints on the shape function 
S), which are independent of the leptonic variable P x . The fact that the total decay rate T u 
is proportional to mjj is not in contradiction with (23). It is instructive to demonstrate how 
these five powers of m b are recovered in our approach. At tree level and leading power the 
functions 3 r 2 and 5F 3 vanish, while 3"i = S{P + ). Integrating over the full range of Pi and P_ 
builds up five powers of (M B — P + ). For the purpose of illustration, let us rename the P + 
variable to Cj in the last integration, so that the total decay rate is given as 



r<1 It/ 12 rM B (~i2 It/ 12 pM B -m b 

r u = ^f/ Q du{M B -ufS{u) = ( ^^ J ^ du{m b + ufS{u) 



2 IT/ . 12 



Gi\V. 



F\ v ub\ / / ,\\5 



1927T 3 



(m b + (lo)Y 



1 + 



1 



(27) 



At tree level, the first moment of the shape function S(u) vanishes. Beyond tree level this is 
no longer the case, and the average (cu) depends on the size of the integration domain. The 
above observation motivates the use of the shape-function scheme [9], in which the 6-quark 
mass is defined as mf F = m£ olc + (u) + 0(l/m b ). After this is done, (27) recovers the form of 
the conventional OPE result. 

Eq. (23) and the above argument tell us that the differential rate is a priori rather insensitive 
to the 6-quark mass in the endpoint region, where P + (and therefore {Cj)) is a small quantity. 
Only when the rates are integrated over a sufficiently wide domain, so that shape-function 
integrals can be approximated using a moment expansion, a dependence on m b enters indirectly 
via the first moment of the leading-order shape function. Likewise, a dependence on other 
HQET parameters such as Ai enters via the sensitivity to higher moments. 

In the remainder of this section we present the various contributions to the structure 
functions 5Fj, following the same line of presentation as we did in the case of B — > X s ^ 
decay in Section 2. As before, while the resulting expressions are "optimized" for the shape- 
function region, they can be used over the entire phase space and give the correct result for 
the total decay rate up to corrections of 0(a 2 ). In the shape-function region, where P + is a 
small quantity, one may organize each 3^ as a series in inverse powers of l/(Mg — P + ). No 
assumption about the variable y is made, which is treated as an 0(1) quantity. 2 A preview 
of the results of our calculation is depicted in Figure 1, which shows an illustration of our 
prediction for the distribution of events in the plane (P + ,P_). 



3.1 Leading-power factorization formula 

The leading-power expressions for the hadronic structure functions Wi have been calculated in 
[9] at one-loop order in renormalization-group improved perturbation theory. At this level 

2 In the shape- function region, where P + <C P_, we have y w p_/m b , which is the variable used in the 
leading-power analysis in [9]. 
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>_ [GeV] 



Figure 1: The hadronic phase space in P + and P_. The light gray region contains background 
from B — > X c l~v decays, while the dark gray region is only populated by B — > X u l~v events. 
The line separating the two regions is the contour where M\ = P+P- = M^. Each point 
represents a B — > X u l~v event in a Monte-Carlo simulation using the results of this paper. 
While the shape-function region of large P_ and small P + is highly populated, there is not a 
single event with P + larger than 3 GeV out of the 1300 events generated. 



does not obtain a contribution, whereas 9^ and 5F3 do. Symbolically, they take the factorized 
form H ui J <S> S, consisting of hard functions H ui and the convolution of the jet function J 
with the leading shape function S. More precisely, 

f p + 

3f\p+,y) = H ui (y,fj, h ) dujym b J(ym b (P + - u), fi { ) S(u, fa) , (28) 

Jo 

where the hard functions are given by 



H ul (y,n h ) = 1 -+ 



C F a s (jji 



h 



47T 



2 



, oi/mj , ymh , 21ny , . it 

-41n 2 ^ + 10 In — - 41ny -4L 2 (l-y) 12 

1-y 6 



H u3 {y, Hh) = , (29) 

47r 1 — y 

and H U 2 = 0. As before, the differential decay rate is independent of the matching scales 
Hh ~ m b an d fa ~ ^ym b AQ CI ). The jet function J has already been given in (7). Note that the 
6-quark mass appears only as the argument of logarithms, where it plays the role of setting 
the renormalization scale. 
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3.2 Kinematical power corrections 



As in the case of B — > X s ^ decay, there is a class of power corrections to the B — > X u l~u decay 
distributions which are small only because of the restriction to certain regions in phase space, 
but which are not associated with new hadronic parameters. In the present case, these terms 
can be extracted from the one-loop expressions derived in [22]. They are then convoluted with 
the leading shape function. As previously, the scale separation that can be achieved for these 
power-suppressed terms is only approximate, and we thus assign a coupling a s (/x) with them, 
where the scale p, is expected to be of order Hi ~ ^m^Qcxy. 

The resulting expressions for the structure functions can be written in a compact form in 
terms of the variables x and y defined in (9) and (25). We find 



C F a s (fx) 



x 



P+ 4tt 
fi(x,y) 



f 

Jo 



dCo S(cu, fii) 

2gi(x,y) 



(1 + x) 2 y(x + y) x(l + x) 2 y 2 (x + y) 
1 f, '>,.(//) 



\ x/ x V xJ 



M B - P + 4tt 

f2(x,y) 



x 



du S(lu, Hi) 
2x g 2 (x,y) 



grkin 



(1 + x) 2 y 2 (x + y) (1 + x) 2 y 3 (x + y) 
1 C F a s (fl) r? + 



In 



(> + !) 



x 



P + 4vr 
fs(x,y) 



+ 



du S(u, Hi) 
2gs(x,y) 



[1 + x) 2 y 3 (x + y) (1 + x) 2 y A (x + y) 
where the functions f iy g; L are given by 



In 1 + ^ 



(30) 



h(x,y) 
9i(x,y) 
f2(x,y) 
92{x,y) 
h{x,y) 



-9y + I0y 2 + x(-16 + 12y + 6y 2 ) + x 2 (l3y - 12) 



-2y 3 - 2xy 2 (4 + y) - x 2 y(12 + Ay + y 2 ) - 4x 3 (y + 2) + 3x 4 (y - 2) , 
y 2 + xy(S + Ay + y 2 ) + 3x 2 y(10 + y) + x 3 (l2 + 19y) + 10x 4 , 
2y 2 + Axyil + 2y) + x 2 y(18 + by) + 6x 3 (l + 2y) + 5x 4 , 
2y 3 (2y - 11) + xy 2 (-94 + 2% + 2y 2 ) + 2x 2 ?/(-72 + 18y + I3y 2 ) 
+ x 3 (-72 - 42y + 70y 2 - 3y 3 ) - 10x 4 (6 - 6y + y 2 ) , 
g 3 (x, y) = Ay 4 - 6x(y - 5)y 3 - Ax 2 y 2 (-20 + 6y + y 2 ) + x 3 y(90 - lOy - 28y 2 + y 3 ) 

+ x 4 (36 + 36y - 50y 2 + Ay 3 ) + 5x 5 (6 - 6y + y 2 ) . (31) 

The above formulae are the exact 0(a s ) corrections to the leading-power expression. This 
means that, when integrated over the entire phase space, they will give rise to the correct 
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result for the total rate up to that order. In the shape-function region (where F + < P_) the 
integrands in (30) can be expanded in powers of 1/mj, by counting y = 0(1) and x = 0(l/ra&). 
Note that this organizes the l/m& expansion as an expansion in powers of the hadronic variable 
1/(Mb — P+). The leading terms read 



?? n{1) (P + ,y) 



J kin(l) | 



p+,v) 



r kin(l) 



(P + ,y) 



1 




C F a s (fi) 


M B - 


P+ 


An 


1 




C F a 8 {p) 


Mb — 


P + 


An 


1 




C F a 8 {p) 


M B - 


P + 


An 



f p+ « 

/ dw S(u, Hi) 
Jo 

°+ 

dw S(w, Hi) 

f P+ ~ 
/ dw S(w, Hi) 

Jo 



6-« + (H_4W 

y \y J x . 



\ 22 8 i y 

4 H — la — 



(32) 



y y x 

Further accuracy can be achieved by adding the next-order corrections, for which we obtain 



?f n(2) (P + ,y) 



1 C F a s (jj) 



[M B - P + f An 



f 

Jo 



dw (P + — lj) S(w, Hi) 



x 



16 3 /12 20 \ , y 
■12 + — + — + — + 61n^ 

y y \y y J x_ 



^ iQi2) (P + ,y) = 



C F a s (fi) 



(M B - P + ) 2 An Jo 



f P + 

/ dCo (P 4 
Jo 



Co) S(w,fj>i 



x 



?3 iQ(2) (P + ,y) 



2 7 A v 
l + - + ^- ^ln- 

y y y x 



C F a s (n) 



(M B - P+) 2 An 



duo (P + — Co) S(w, Hi) 



o 



x 



69 64 /52 28\ , y 
-6 + + — ln- 

y yl y J 3. 



(33) 



In the various phase-space regions of interest to the determination of \V u b\, the above terms 
(32) and (33) approximate the full result (30) very well (see Section 5 below). 

Let us comment here on a technical point already mentioned in the Introduction. When 
combined with the phase-space factors in (23), the exact expressions for S^ 111 in (30) are regular 
in the limit P — > P + , corresponding to y — > 0. However, this feature is not automatically 
ensured when the structure functions, but not the phase-space factors, are expanded about the 
heavy-quark limit. With our choice of the variables x and y, we encounter terms as singular 
as l/y n at n-th order in the expansion, as is obvious from the explicit expressions above. 
Phase space scales like y 2 in the limit y — > (note that Pi — > P + as P_ — > P + because of 
(2)), so that the results (32) and (33) can be applied without encountering any kinematical 
singularities. In order to achieve this, it was crucial to define the variable y in the way we 
did in (25). We emphasize this point because straightforward application of the technology of 
SCET and HQET developed in [16, 31, 32] would give an expansion of the structure functions 
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Hi in powers of whereas phase space is proportional to 4p 2 = (p_ — p + ) 2 oc y 2 . In the 

kinematical region where p + < 0, which is allowed due to off-shell effects in the B meson, this 
leads to singularities as p_ — > 0. In order to avoid these singularities, we have reorganized 
the SCET expansion as an expansion in — p+) instead of l/p~, where \p + \ <C p_ in the 

shape-function region. 



3.3 Subleading shape- function contributions 



The contributions from subleading shape functions to arbitrary B — > X u Pu decay distri- 
butions have been derived (at tree level) in [16, 31, 32]. The results involve the same set 
of subleading shape functions as previously discussed in Section 2.3. Again, the structure 
function 5F 2 does not obtain a contribution, while 

u(P + )-v(P + 



l 



(A-P + )S(P + )+t(P + ) + 



J: 



hadr(l) 



p+,v) 



2 



-(A-P + )S(P + )-2t(P + ) + 



i(p + ) + v(p + ) 



y 



(34) 



At this point we recall the discussion of Section 2.3, where we have argued that the B — > X s 7 
photon spectrum should be used to fit the function S of (14), which is defined to be a linear 
combination of the leading shape function S and the subleading shape functions t, u, v. 
When the above results are rewritten in terms of the new function S nothing changes in the 
expressions for 9 r ^°' ) except for the simple replacement S — > S, which we from now on assume. 
At the level of subleading shape functions 3 r £ adr ^ = and g^adr(i) remain unchanged, while 



^ dr(1) (P + ,y) = 



Mn-P+ 



-(A - P + ) S(P+) + 2t(P + ) + [u(P + ) - v(P + )\ - 



(35) 



It follows that there reside some linear combinations of subleading shape functions in the triple 
differential decay rate that cannot be extracted from information on the photon spectrum in 
B — > X s 7 decays. In the end, this dependence gives rise to a theoretical uncertainty. 



3.4 Residual hadronic power corrections 

In analogy with our treatment for the case of B — > X s ^ decay, we start from the expressions 
for the I /m 2 , corrections to the triple differential B — > X u l~v decay rate obtained by applying 
the OPE to the hadronic tensor [12, 13]. Converting these results into the (v,n) basis and 
changing variables from v ■ q and q 2 to p + = n ■ p and p- = n ■ p, we find 

.1,(2) „ J, , 2X 1 -3X 2 \ /2AJ-3A2 5Ai + 15X 2 \ ^ 

= Hp + ) (i + ) + « (p + ) {—_ ^— )-'W T . 

W? = S (P + ) ("^^) • (36) 

1 w? + WP + i Mf » = *±> ( 2Xl + 12Aa - 4A ' + 9A ^ + f^tl I ^ + 4A 2 ) 
4 j/ P- \ 3P- 3m 6 y V 3 / 
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The desired 1/(M B — P + ) 2 corrections to the structure functions 3^ can then be extracted 
by expanding the leading and subleading contributions 3^ and 5F'! iadr ' 1 ) j n terms of their 
moments in (18), and by subtracting the results from (36). Following the same procedure 
as in Section 2.4 to express the remaining power corrections in terms of the leading shape 
function, we obtain 

3->»dr(2), p , 1 / 4A 1+ 24A 2 4A!+9A 2 \g 

^ (Af B -P + )4 3y* 3y ) S(P+h (3?) 

These expressions remain unchanged when the shape function S is used instead of S. 



3.5 Weak annihilation contributions 

In the OPE calculation several contributions appear at third order in the power expansion: 
I/mi, corrections to the kinetic and chromo-magnetic operators, the Darwin and spin-orbit 
terms, and weak annihilation contributions. The Darwin and spin-orbit terms correspond to 
the forward 5-meson matrix elements of (light) flavor- singlet operators [40]. The correspond- 
ing HQET parameters p 3 D and p\ s can in principle be extracted from moments of inclusive 
B — > X c l~u decay spectra. They are insensitive to the flavor of the spectator quark inside 
the B meson. The weak annihilation contribution, on the other hand, results from four-quark 
operators with flavor non-singlet structure. Graphically, this contribution corresponds to a 
process in which the b and u quark annihilate into a W~ . Weak annihilation terms come with 
a phase-space enhancement factor of 16tc 2 and so are potentially more important than other 
power corrections of order l/mf. Because of the flavor dependence, these contributions can 
effect neutral and charged B mesons differently [41]. One choice of basis for the corresponding 
four-quark operators is [42] 

f 2 M 2 - - f 2 /If 2 

(B\b Ll »u L u Ll »b L \B) = } JL_tL Bl , (B\b R u L u L b R \B) = J -^- B 2 , (38) 

where fs is the B- meson decay constant, and Bi are hadronic parameters. In the vacuum 
saturation approximation they are given by B\ — B% — 1 for charged B mesons and B\ — B% — 
for neutral ones. The total semileptonic rate is proportional to the difference (B 2 —Bi), which 
implies that the weak annihilation contribution would vanish in this approximation. Currently, 
only rough estimates are available for the magnitude of the deviation of this difference from 
zero. The resulting effect on the total branching ratio is [43] 

SB(B X., I-,) » 3.9 (^) 2 (^i) |VU» . (39) 
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Again, we expect this effect to be different for charged and neutral B mesons. The most 
important feature of weak annihilation is that it is formally concentrated at the kinematical 
point where all the momentum of the heavy quark is transferred to the lepton pair [41]. At the 
parton level this implies that the corresponding contribution is proportional to S(q 2 — rn 2 ). It 
is therefore included in every cut that includes the q 2 endpoint, and its effect is independent 
of the specific form of the cut. 

We suggest two different strategies to control this effect. The first is to include it in the error 
estimate as a constant contribution proportional to the total rate. A recent study [44] puts a 
limit on this effect of ±1.8% on the total rate (at 68% confidence level) by analyzing CLEO 
data. The second one is to impose a cut q 2 < q 2 nax , thus avoiding the region where the weak 
annihilation contribution is concentrated. The maximal value of q 2 is (M B — M^) 2 , but one 
must exclude a larger region of phase space, such that the excluded contribution to the decay 
rate at large q 2 (corresponding to a region near the origin in the (P_, P + ) plane) can be reliably 
calculated. In our numerical analysis, we will study the effect of a cut q 2 < (Mb — Mp) 2 , 
which satisfies this criterion. 

For completeness, we note that even after the weak annihilation contribution near maxi- 
mum q 2 has been removed, there could in principle exist other, flavor-specific contributions to 
the semileptonic decay amplitudes that are different for charged and neutral B mesons. The 
leading terms of this kind contribute at order 1/rrib in the shape-function region and are pa- 
rameterized by a set of four-quark subleading shape functions [16, 31, 32]. Model estimates of 
these contributions show that they are very small for all observables considered for an extrac- 
tion of \ V u b\ [32, 45]. If only flavor- averaged decay rates are measured, the effects of four-quark 
subleading shape functions can be absorbed entirely by a redefinition of the functions u(Co) 
and v{Cj) [16], without affecting the moment relations in (18). 

4 Shape-function parameterizations 

Hadronic-physics effects enter the description of inclusive decay rates via non-perturbative 
shape functions. Perturbation theory cannot tell us much about the local form of these 
functions, but moments of them are calculable provided that the domain of integration is 
much larger than Aqcd- Since the shape functions contain information about the internal 
structure of the B meson, knowledge of them relates directly to the determination of the 
6-quark mass m b , the kinetic-energy parameter Ai, and in principle the matrix elements of 
higher- dimensional operators. Improved measurements of the shape of the B — > X s 7 photon 
spectrum will therefore lead directly to a more precise determination of HQET parameters. 
This argument can be turned around to constrain the leading shape function using knowledge 
of m b and Ai from other physical processes such as a b — > c moment analysis [46]. We 
emphasize, however, that there are obviously infinitely many locally different functions that 
have identical first few moments. In this section we present a few functional forms that can 
be used to model the shape functions and to fit the current experimental data. 

To achieve stringent constraints on the leading shape function a precise definition of the 
HQET parameters is required. It is a well-known fact that the pole-mass scheme introduces 
uncontrollable ambiguities. To avoid these uncertainties several short- distance definitions 
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have been proposed, such as the MS scheme, the potential-subtraction scheme [47], the T(IS') 
scheme [48], the kinetic scheme [49], or the shape-function scheme [9]. While the decay rates 
are of course independent of the particular choice, it is advantageous to use a mass scheme 
that is designed for the physics problem at hand. In the case of inclusive B decays into light 
particles, this is the shape-function scheme. 



4.1 Models for the leading shape function 

Model-independent constraints on the shape function S(lu, Hi) can be derived by analyzing 
moments defined with an upper limit of integration cu , i.e. 



M N (u ,fii) 



rui 



(40) 



For practical applications, ojq should be taken of order the size of the window where the B — > 
Xg'-f photon spectrum is experimentally accessible, Cjq = Mb — 2E" un with E™ m 1.8 GeV. 
These moments can be expanded in terms of matrix elements of local operators as long as ujq 
is large compared to Aqcd- In the shape-function scheme, HQET parameters are defined to 
all orders in perturbation theory through ratios of such moments, e.g. [9] 



M 1 (fj,f + A(fj,f,fj, i ),fi i ) 
M (/i/ + A^/,^),^) 

M 2 {fif + A{fif,fii),fii) 
M (nf + A(n f ,Hi),Hi) 



(41) 



Here, the factorization scale /// >> Aqcd is related to the size of the integration domain via 
the implicit equation ujq = fif + A(/ij,/ij). In practice Hf is close to the intermediate scale 
Hi- At tree level, the relations between parameters in the shape- function scheme and the pole 
scheme are A(Hf,Hi) = A polc and H^(Hf,Hi) = — Ai. The corresponding relations at one- and 
two-loop order have been worked out in [9] and [50] , respectively. These relations allow us to 
obtain precise determinations of A(/j/,/ij) and /4(/t/,/tj) from other physical processes. 

For reference purposes, it is helpful to quote values for A and h\ using only a single scale 
/j* instead of two independent scales /i/ and /ij. To one-loop order, these parameters can be 
related to those determined from the moments via [9] 



7/ \ 7/ \ | C F a s (n*) C F a s {ni) 



71 



7T 



1-2 1 



1 C F a s {^) _ C F a s (fn) f 1 3 ^ 

2ir it \2 Hi 



3/4 
+ 3^ 



ln^ 

Hi 



(42) 



2 C F aJjM) ln Hi 

7T Hi 



where we have neglected higher-dimensional operator matrix elements that are suppressed 
by inverse powers of /t/. A typical choice for the scale h* is 1-5 GeV, which we will use as 
the reference scale throughout this work. It will be convenient to connect the parameter A 



19 



extracted from the first moment of the shape function with a low-scale subtracted quark-mass 
definition referred to as the "shape- function" mass. Following [9], we define 

™>bbJ>f,(J>i) = M B -A(nf,Hi). (43) 

The general procedure for modeling the leading shape function S(ui, Hi) from a given func- 
tional form F(u) is as follows. The shape of F(cu) is assumed to be tunable so that it can be 
used to fit the B — > X s 7 photon spectrum. Only the norm of the shape function is fixed the- 
oretically. Note that the moment relations (41) are insensitive to the norm, so that formulae 
for A and h\ follow directly from the functional form of F(u). Examples of such formulae will 
be given below. We define moments m\P(u q ) of F in analogy with (40). The first relation in 
(41) implies that for a given u the factorization scale is 

M^Vo) 

Now that /if is known, the norm is determined by requiring that the zeroth moment of the 
shape function is [9] 

(45) 

IF] " ~ 

It follows that [M q (lj , fi^/M^ (u>o)] F{Q) serves as a model of S(cu,Hi) or S(uj,Hi)- 

We now suggest three two-parameter models for the leading-order shape function based 
on an exponential-type function F^- Cxp \ a gaussian-type function _p(s auss ) ; and hyperbolic-type 
function F^ hyp ^. We use two parameters that can be tuned to fit the photon spectrum: a 
dimensionful quantity A which coincides with the position of the average (a)), and a positive 
number b which governs the behavior for small Co. The functions we propose are 

F <^ ;A>6H ^g) 6 ^_^ p) £), 



7r \ pLi Hi 24 J it \ Hi 2 J 3/j, f 



F<*°-\u,; A, b) = (£) exp (-rf (gauss) ^) , 

F^(u;Kb) = ^ 'cosh- 1 (d(hyp)|) • (46) 

For convenience, we normalize these functions to unity. The parameters are determined 
by the choice A = (a)). We find 

^ = r(6) ' ^(cxp) — b , 

2d b/2 /r(±±b-\\ 2 

^(gauss) = ^s) j , LL2J > (4?) 



7\r(hyp) 




r(6/2) ' 

[4^ (hyp) ] fc 6 C(l + 6^)-C(l + &, j) 

2F(6) [C(6, 1) - C(6, !)] ' ^ 4 C(6, J) - C(6, !) 
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Figure 2: LEFT: Different functional forms for the leading shape function. We show 
F^(u, A,2) (solid), Ffe auss )(^,A,2) (dotted), and F (hyp) (a), A, 2) (dash-dotted) as functions 
of the ratio u/A. Right: The same functions with the parameters A and b tuned such that 
m&(/i*, //*) = 4.61 GeV and /i 2 (//*,//*) = 0.2 GeV 2 . See text for explanation. 



where ((b,a) = YlkLoi^ + a ) ^ s ^ ne generalized Riemann zeta function. An illustration of 
the different functional forms is given on the left-hand side in Figure 2. We show a plot with 
the choice 6 = 2, corresponding to a linear onset for small u. 

For the first two models, analytic expressions for the HQET parameters A and /i 2 are 
available. Following the discussion above, we compute the moments on the interval [0, ci>o] and 
find for the exponential form F^ exp \u; A, b) 



A(fi f ,fii 

where fif = u — A(/xj, fi/j 
A(/i/,/ij) = 



a r(i + &) 
b 



r(i + b, b -f) 



r(b, m 



m ^, A 

a 2 r(2 + b) - r(2 + b, 
1? r(6) -r(6,^|a; 



1^0 ' 

A . 



A((J, f ,fJ,i) 



(48) 



A similar calculation for the gaussian form F^ auss \u; A, b) yields 
A 



I-!/ 1+6 \ T-i/l + b ^(gauss)^O ' 

1 y.—) - 1 \—-> — — . 



^/^(gauss) 

A 2 



r(|: 



r(i 



A 2 



rri 



b ^(gauss)^o - 
2' A 2 



(gauss) 



r(I) 



p/ b "(gauss) ' 

1 ^ 2 ' A 2 



(49) 



The corresponding relations for F( hyp )(cj; A, 6) must be obtained numerically 

Ultimately the shape function should be fitted to the B — > X s -y photon spectrum, and 
the above equations then determine A and /x 2 . On the other hand, these formulae can be 
inverted to determine A and b from the current values of the HQET parameters. For example, 
if we adopt the values mf,(//*, = 4.61 GeV and /x 2 (//*, /i*) = 0.20 GeV 2 for the parameters 
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in (42) at //* = 1.5 GeV, then we find the parameter pair A ss 0.72 GeV, b ps 3.95 for the 
exponential model, A 0.71 GeV, b ~ 2.36 for the gaussian model, and A pa 0.73 GeV, b ps 3.81 
for the hyperbolic model. On the right-hand side of Figure 2 we show these three different 
functions plotted on the interval [0, ujq] over which the moment constraints are imposed. While 
the exponential (solid) and hyperbolic (dash-dotted) curves are barely distinguishable, the 
gaussian model has quite different characteristics. It is broader, steeper at the onset, faster to 
fall off, and the maximum is shifted toward larger Co. 

An important comment is that, once a two-parameter ansatz is employed, the shape- 
function parameters (i.e., m b and //^) can either be determined from a fit to the entire photon 
spectrum, or to the first two moments of the spectrum. Both methods are equivalent and 
should yield consistent results. If they do not, it would be necessary to refine the ansatz for 
the functional form of the shape function. 

In most applications shape functions are needed for arguments Co of order Aqcd- However, 
in some cases, like the ideal cut on hadronic invariant mass, Co is required to be as large as 
Md, which is much larger than Aqcd- The large-u) behavior of the shape functions can be 
computed in a model-independent way using short- distance methods. For the leading shape 
function, one finds [9] 

$(£, » Aqcd, m ) = - "^ . 1 , U k ±± + l) + . . . . (50) 

71 UO - A V A*Z / 

Note that this radiative tail is negative, implying that the shape function must go through 
zero somewhere near Co ~ few Aqcd- For practical purposes, we "glue" the above expression 
onto models of the non-perturbative shape function starting at Co = A + Hi/y/e ~ 1.6 GeV, 
where the tail piece vanishes. In this way we obtain a continuous shape-function model with 
the correct asymptotic behavior. We stress that for applications with a maximal P + not larger 
than about 1.6 GeV the radiative tail of the shape function is never required. This includes 
all methods for extracting \V u b\ discussed later in this work, except for the case of a cut on 
hadronic invariant mass, M x < M , if M is above 1.6 GeV. 

4.2 Models for subleading shape functions 

In the last section we have been guided by the fact that the B — > X s 7 photon spectrum is 
at leading power directly determined by the leading shape function. This helped in finding 
models that have roughly the same shape as the photon spectrum. At the subleading level 
considered here, however, no such guidance is provided to us. The available information 
is limited to the tree-level moment relations (18), stating that the norms of the subleading 
shape functions vanish while their first moments do not. In [16], two classes of models have 
been proposed, in which the subleading shape functions are "derived" from the leading shape 
function. A particularly simple choice is 

2A 

i(Co) = -X 2 S'(Co), u(Co) = —^S'(Co), v(Co) = X 2 S'(Co). (51) 

o 

Below, we will sometimes refer to this set of functions as the "default choice" . We choose the 
parameter — Ai in the expression for u{Co) (as well as in the expressions for the second-order 



22 



hadronic power corrections) to coincide with the quantity fi 2 (fif,fii) given in (48) and (49). 
However, for consistency with the tree-level moment relations, we identity the parameter A in 
(15) and (35) with the quantity A(/ij,/ij) evaluated in the limit where o;o — > oo. This implies 
A = A for all three types of functions and ensures that the subleading shape functions have 
zero norm when integrated over < u < oo. 

There are of course infinitely many possibilities to find models for subleading shape func- 
tions that are in accordance with (18). Any function with vanishing norm and first moment 
can be arbitrarily added to any model for a subleading shape function without violating the 
moment relations. Several such functions have been proposed in recent work on subleading 
shape functions, see e.g. [16, 30, 32, 45]. Specifically, we define the functions 
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(52) 



where z = Cj/VLq, and the reference quantity VLq = O(Aqcd) depends on the type of function, 
namely Q = A for hi and h 2 , — | A for h 3 , and fi = ^ A for h^. The quantity a is a 
free parameter. The functions (52) have by construction vanishing norm and first moment. 
Their second moments are given by the parameter M 2 , provided the normalization constants 
are chosen as = 1 for hi and h 2 , and 



N =1- 



71 



2(tt-2) a ' 



N = 1- 



7T 2 (4-7T) 1 

8(71- -2) a 



(53) 



for h 3 and h±, respectively. The values for the parameters a and M 2 should be chosen such 
that the following characteristics of subleading shape functions are respected: First, they are 
dimensionless functions, so that their values are naturally of 0(1) for Co ~ Aq CD . Secondly, 
when integrated over a sufficiently large domain, their contributions are determined in terms of 
their first few moments. In particular, this implies that for values of Co ^> Aqcd the integrals 
over the subleading shape functions must approach zero. Taking these considerations into 
account, we use M 2 = (0.3 GeV) 3 in all cases and choose a = 3.5 for hi, a — 5 for h 2 , and 
a = 10 for h 3 and /14. 

Given the four functions (52), we can construct several new models for the subleading 
shape functions t(u), u{Co), and v(Co). For each function, we construct a set of 9 models by 
adding or subtracting any of the functions h n (uj) to the default choice in (51). Together, this 
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Figure 3: Nine models for the subleading shape function u{Co) obtained by adding or subtract- 
ing one of the four functions h n (Co) to the default model in (51), shown as a thick line. See 
text for explanation. 



method yields 9 3 = 729 different sets {tj(u>), Uj(co), with i,j,k = 1,...,9. This large 

collection of functions will be used to estimate the hadronic uncertainties in our predictions 
for partial decay rates. Note that for most of these sets we no longer have U(Co) = —Vk(Co), 
which was an "accidental" feature of the default model (51). The fact that the two functions 
have equal (but opposite in sign) first moments does not imply that their higher moments 
should also be related to each other. 

For the case of u(Co) the resulting functions are shown in Figure 3, where we have used 
the exponential model (46) with parameters A = 0.72 GeV and b = 3.95 for the leading shape 
function. In the region Co ~ Aqcd they differ dramatically from each other, while the large Co 
dependence is dominated by the moment relations (18). 

4.3 Illustrative studies 

We stressed several times that the calculation of the hadronic tensor is "optimized" for the 
shape-function region of large P_ and small P + , while it can smoothly be extended over the 
entire phase space. The notions "large P_" and "small P + " are to be understood as the sizes 
of integration domains for P_ and P + . Only when the differential distributions are integrated 
over a sufficiently large region in phase space, global quark-hadron duality ensures that the 
partonic description used in the present work matches the true, hadronic distributions with 
good accuracy. A more ambitious goal would be to calculate the differential decay rate point 
by point in the (P + , P_) plane. This can be done invoking local quark-hadron duality, as long 
as there is a sufficiently large number of hadronic final states contributing to the rate at any 
given point in phase space. 

It is instructive to integrate the triple differential decay rate (23) over the leptonic variable 
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Figure 4: Left: Theoretical prediction for the double differential decay rate. The light area 
represents a large decay rate. Black regions denote areas where the decay rate is close to zero. 
The dotted line is given by P+P_ = Ml, which means that charm background is located in the 
upper wedge. See text for further explanation. Right: The P + spectrum extended to large 
values of P + . The thin solid line denotes the leading-power prediction, the dashed line depicts 
first-order power corrections, the dash-dotted line shows second-order power corrections, and 
the thick solid line is their sum. 

Pi in the range P + < Pi < P_, which yields the exact formula 

P + f 

P_)J 2 + (P_-P + )? 3 }. (54) 

Our theoretical prediction for the double differential decay rate (54) is shown on the left-hand 
side of Figure 4. We use the exponential model for the leading shape function with parameters 
m&(//*, f/,*) = 4.61 GeV and /U*) = 0.2 GeV 2 , as well as the default choice (51) for the 

subleading shape functions. For very small P_ values the rate turns negative (to the left of the 
gray line in the figure), signaling a breakdown of quark-hadron duality. It is reassuring that 
the only region where this happens is the "resonance region", where the hadronic invariant 
mass is of order Aq CD , and local duality breaks down. 

Another useful quantity to consider is the differential P + rate, which is obtained by inte- 
grating the double differential rate over P_ in the range P + < P_ < M B - The resulting P + 
spectrum is shown on the right-hand side of Figure 4. In the plot we also disentangle the 
contributions from different orders in power counting. 
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5 Predictions and error estimates for partial rates 



Before discussing predictions for partial B — > X u l"v rates for various kinematical cuts, let 
us recapitulate the ingredients of the calculation and general procedure. We have presented 
expressions for the triple differential decay rate, which can be organized in an expansion in 
inverse powers of {Mb — P+)- The leading-power contribution is given at next-to- leading 
order in renormalizat ion-group improved perturbation theory. At first subleading power two 
contributions arise. The first type involves subleading shape functions and is included at tree 
level, while the second type contributes perturbative corrections of order a s that come with 
the leading shape function. Further contributions enter at second subleading power and are 
again of the two types: perturbative corrections of order a s and non-perturbative structures 
at tree level. In summary, then, partial rates can be computed term by term in an expansion 
of the form 

r u = rf) + (rfw + r* adr «) + (rf 1 ^ + rj; adr ( 2 )) + . . . . (55) 

The goal of this section is to test the convergence of this series expansion and to perform 
a thorough analysis of uncertainties. For the kinematical corrections T^ n ^ the sum of all 
terms is known and given by the expressions in (30), while the first two terms in the series 
correspond to the expanded results in (32) and (33). We will find that in all cases of interest 
the first two terms give an excellent approximation to the exact result for T^ 111 . 

For the purpose of illustration, we adopt the exponential model for the shape function and 
present numerical results for two sets of input parameters, which are biased by the results 
deduced from fits to B — > X c l~v moments [50]. Specifically, we use m b (//*, //*) = 4.61 GeV, 
= 0.2 GeV 2 (set 1) and m b (//*,//*) = 4.55 GeV, /i 2 (/i*, //*) = 0.3 GeV 2 (set 2). The 
values of the 6-quark mass coincide with those obtained at two-loop and one-loop order in 
[50] (see also the discussion below), while the values of /i 2 are close to the corresponding 
values in that reference. As was mentioned before, in the future the leading shape function 
S{Cj,jj,i) should be extracted from a fit to the B — > X s ^ photon spectrum, in which case the 
uncertainty in its shape becomes an experimental error, which can be systematically reduced 
with improved data. In the process, the "theoretically preferred" parameter values used in 
the present work will be replaced with the "true" values extracted directly from data. While 
this will change the central values for the partial rates, our estimates of the theoretical errors 
will only be affected marginally. 

The different sources of theoretical uncertainties are as follows: First, there are uncertain- 
ties associated with the functional forms of the subleading shape functions. To estimate them, 
we take the spread of results obtained when using the large set of different models described 
in Section 4.2, while the central value for a partial decay rate corresponds to the default model 
(51). Secondly, there are perturbative uncertainties associated with the choice of the match- 
ing scales /i/j, /ij, and p. Decay rates are formally independent of these scales, but a residual 
dependence remains because of the truncation of the perturbative series. Our error analysis 
is as follows: 

• The hard scale fih is of order m b . In perturbative logarithms the scale appears in the 
combination {ym^/ fJ-h), see e.g. (29). To set a central value for fih we are guided by the 
average (y)rrib. The leading term for the double differential decay rate d 2 T u /dP + dy is 
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proportional to 2y 2 (3 — 2y). It follows that the average y on the interval [0, 1] is 0.7. 
However, in some applications y is not integrated over the full domain. Also, there are 
large negative constants in the matching correction if ul in (29), whose effect can be 
ameliorated by lowering the scale further. In the error analysis we use the central value 
of ii h = m b /2 2.3 GeV and vary the scale by a factor between 1/a/2 and For the 
central value a s (fx h ) « 0.286. 

• The intermediate scale fa ~ a/^Aqcd serves as the renormalization point for the jet 
and shape functions. We fix this scale to fa = 1.5 GeV. Variations of fa would affect 
both the normalization and the functional form of the shape function, as determined 
by the solution to the renormalizat ion-group equation for the shape function discussed 
in [9, 11]. In practice, effects on the shape are irrelevant because the shape function 
is fitted to data. The only place where the intermediate scale has a direct impact on 
the extraction of |T4b| is through the normalization of the shape function (45). In the 
analysis we therefore estimate the uncertainty by assigning the value ±( " s ^ i ) ) 2 as a 
relative error, where a s (fa) ~ 0.354. 

• The scale p, appears as the argument of a s in the perturbative contributions r^ m . We 
vary p, from fa/ \f2 to \[2fa with the central value p, — fa — 1.5 GeV. 

These three errors are added in quadrature and assigned as the total perturbative uncertainty. 
Finally, we need to estimate the effects from higher-dimensional operators at third and higher- 
order in power counting. If the considered cut includes the region of phase space near the 
origin (P + ~ P_ ~ Aqcd)j then the dominant such contributions are weak annihilation effects, 
which we have discussed in Section 3.5. From the analysis in [44] one can derive a bound on 
the weak annihilation contribution that is ±1.8% of the total decay rate, for which we take 
T„ ~ 70 iKftPps -1 (see below). The resulting uncertainty 5T™ A = ±1.3 |Kb| 2 ps _1 affects all 
partial rates which include the region near the origin in the (P+, P_) plane. The uncertainty 
from weak annihilation can be avoided by imposing a cut q 2 < g^ ax (see Section 5.6). For 
all observables considered in the present work, other power corrections of order l/m\ can be 
safely neglected. This can be seen by multiplying the contributions from second-order hadronic 
power corrections to the various decay rates (called rjj adr ( 2 )) D y an additional suppression factor 
Aqcd/w 6 ~ 0.1. 

The following subsection contains a discussion of the total decay rate. In the remainder of 
this section we then present predictions for a variety of kinematical cuts designed to eliminate 
(or reduce) the charm background. These partial rates can be computed either numerically or, 
in many cases, semi-analytically. In Appendix B we discuss how to perform the integrations 
over the kinematical variables p and P_ analytically. 

5.1 Total decay rate 

Before presenting our predictions for the various partial decay rates, it is useful to have an 
expression for the total B — > X u l~v decay rate expressed in terms of the heavy-quark parame- 
ters defined in the shape-function scheme. We start from the exact two-loop expression for the 
total rate derived in [51], add the second-order hadronic power corrections, which are known 
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at tree level [12, 13], and finally convert the parameters m b and Ai from the pole scheme to 
the shape-function scheme. The relevant replacements at two-loop order can be taken from 
[50] and read 



ml° le = m b + 0.424^0^ (jj) 



1 [ 1.357+ 1.326 In — + 0.182^% 1 a s U) 

A** K 



| 3^-^- 0.330^^) | 
2m b 

-Ax =/4 + 0.330/i^ s 2 (/i) + ... , (56) 

where here and from now on mj = m;,^,^) and /x 2 = /i 2 (/i*, /i*) are defined in the shape- 
function scheme. At a reference scale /i* = 1.5 GeV the values of these parameters have 
been determined to be m b = (AM ± 0.08) GeV and /i 2 = (0.15 ± 0.07) GeV 2 [50] , 3 where we 
account for the small l/m b correction to the relation for the pole mass in the above formula 
(corresponding to a shift of about —0.02 GeV in m b ), which was not included in that paper. 
The resulting expression for the total decay rate is 



r. = «^{l + CM (-0.768 + 2.122^) 



+ OL 2 s {n) 

3(/4 " A 2 ) 



(2 \ 2 
1.249 + 2.814 In + 0.386 % ) ^- + 0.811 ^ 
A** Kj rn b ml 



+ ...}. (57) 



We observe that for /i* ~ 1.5 GeV and /i = 0(m b ), the perturbative expansion coefficients 
are strongly reduced compared to their values in the pole scheme (—0.768 and —2.158, re- 
spectively), indicating a vastly improved convergence of the perturbative expansion. For 
m b = fi = 4.61 GeV, and /j, 2 = 0.15 GeV 2 we obtain for the one-loop, two-loop, and power 
corrections inside the brackets in (57): {1 — 0.017 — 0.030 — 0.004}. All of these are very small 
corrections to the leading term. 

Including the uncertainties in the values of m b and /i^ quoted above, and varying the 
renormalization scale \i between m b and m b /2 (with a central value of m b /y/2), we get 

F / m u \ 4.81 



|Kfc| 2 ps 



68.0^| [m b ] T 0.7 [£] ±g;g M = (68.0 =f 0.7 [^] ±g;g M) (^^iv) ' • (58) 



Here and below, we quote values for decay rates in units of \ V ub \ 2 ps _1 . To convert these results 
to partial branching fractions the numbers need to be multiplied by the average S-meson 
lifetime. Without including the two-loop corrections, the central value in the above estimate 
increases to 70.6. For comparison, with the same set of input parameters our new approach 
based on (23) predicts a total decay rate of T u = (71.4 tg; 2 , ± 0.5) | V^,6 1 2 ps _1 , where the first 

3 The values obtained from a one-loop analysis are m b = (4.55 ± 0.08) GeV and [i\ — (0.34 ± 0.07) GeV 2 . 
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Table 1: Partial decay rate T U (E Q ) for a cut on charged- lepton energy Ei > E Q in the £>-meson 
rest frame, given in units of \V u b\ 2 ps _1 . Predictions are based on the shape- function parameter 
values m b = 4.61 GeV, & = 0.2 GeV 2 (top) and m b = 4.55 GeV, /4 = 0.3 GeV 2 (bottom). 



E [GeV] 


Mean 


Subl. SF 


Pert. 


Total 


1.9 
2.0 
2.1 
2.2 
2.3 
2.4 


24.79 

18.92 

13.07 

7.59 

3.12 

0.42 


±0.53 
±0.60 
±0.71 
±0.81 
±0.89 
±1.05 


+1.90 
-1.66 
+1.35 
-1.20 
+0.82 
-0.75 
+0.38 
-0.34 
+0.15 
-0.16 
+0.16 
-0.22 


+2.34 
-2.15 
+1.95 
-1.84 
+ 1.66 
-1.63 
+1.55 
-1.54 
+1.55 
-1.55 
+ 1.65 
-1.65 


1.9 
2.0 
2.1 
2.2 
2.3 
2.4 


21.10 

15.83 

10.73 

6.12 

2.47 

0.29 


±0.53 
±0.60 
±0.68 
±0.74 
±0.84 
±0.99 


+1.57 
-1.35 
+1.08 
-0.94 
+0.64 
-0.55 
+0.31 
-0.23 
+0.17 
-0.22 
+0.18 
-0.24 


+2.08 
-1.92 
+ 1.77 
-1.68 
+1.57 
-1.54 
+1.50 
-1.48 
+1.53 
-1.53 
+1.61 
-1.62 



error accounts for perturbative uncertainties while the second one refers to the modeling of 
subleading shape functions (to which there is essentially no sensitivity at all in the total rate). 
The fact that this is in excellent agreement with the direct calculation using (57) supports the 
notion that the formalism developed in this work can be used to describe arbitrary B — > X u l~v 
decay distributions, both in the shape-function region and in the OPE region of phase space. 

5.2 Cut on charged-lepton energy 

Traditionally, the most common variable to discriminate against the charm background is the 
charged-lepton energy Ei. As long as one requires that E\ is bigger than (M| — Mq)/2Mb ~ 
2.31 GeV, the final hadronic state cannot have an invariant mass larger than M D . For this 
ideal cut, and using the default set of subleading shape functions, we find 

r^ ) + (r!f n{1) ± r!r dr(1) ) + (r!f n{2) ± rJr dr(2) ) 

= [ 6.810 + (0.444 - 3.967) + (0.042 - 0.555)] |K&| 2 ps _1 . 

The corrections from subleading shape functions are quite sizable, in accordance with the 
findings in [28, 29, 30]. Note that the sum 

r kin(l) + r ki„(2) =04g6 

is an excellent approximation 

to the exact result r^ m = 0.482 (all in units of |K(ft| 2 ps _1 ) obtained using (30), indicating that 
the expansion of the kinematical power corrections is converging rapidly. The same will be 
true for all other observables considered below. 
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Table 2: Same as Table 1, but for the partial decay rate 7r„ (Eq) for a cut on lepton energy 
Ei > Eq in the T(4S) rest frame. 



E [GcV] 


Mean 


Subl. SF 


Pert. 


Total 


1.9 
2.0 
2.1 
2.2 
2.3 
2.4 


24.82 

19.00 

13.25 

7.99 

3.83 

1.31 


±0.54 
±0.61 
±0.71 
±0.78 
±0.86 
±0.99 


+1.91 
-1.66 
+1.37 
-1.21 
+0.85 
-0.76 
+0.42 
-0.37 
+0.18 
-0.13 
+0.10 
-0.14 


+2.35 
-2.15 
+ 1.96 
-1.85 
+ 1.68 
-1.63 
+ 1.54 
-1.53 
+ 1.54 
-1.53 
+ 1.61 
-1.61 


1.9 
2.0 
2.1 
2.2 
2.3 
2.4 


21.16 

15.94 

10.94 

6.49 

3.05 

0.98 


±0.54 
±0.60 
±0.68 
±0.74 
±0.84 
±0.92 


+1.58 
-1.35 
+1.10 
-0.95 
+0.66 
-0.57 
+0.34 
-0.26 
+0.17 
-0.18 
+0.13 
-0.18 


+2.09 
-1.93 
+1.78 
-1.69 
+1.58 
-1.54 
+1.50 
-1.48 
+1.53 
-1.53 
+ 1.56 
-1.57 



In practice, the cut on E\ can be relaxed to some extent because the background is well 
understood, thereby increasing the efficiency and reducing the impact of theoretical uncer- 
tainties. Our findings for different values of the cut E are summarized in Table 1. Here and 
below, the columns have the following meaning: "Mean" denotes the prediction for the partial 
decay rate, "Subl. SF" the uncertainty from subleading shape functions, and "Pert." the total 
perturbative uncertainty. In the column "Total" we add the stated errors plus the uncertainty 
from weak annihilation in quadrature. 

Experiments often do not measure the partial rates in the 5-meson rest frame, but in the 
rest frame of the T(AS) resonance produced in e + e" collisions. Boosting to the T(AS) frame 
with P = v/c 0.064 has a small effect on the spectrum and rates. The exact formula for 
this boost is [20] 



7ri T) (£o) 



ft 



i rM B /2 
P- J/3-Eo 



dE 



(3 + — max ( p_ , 



Eq 

~E 



(60) 



where f3± = y/1 ± f3/y/l =F ft and the factor 7 = 1/y 1 — ft 2 ~ 1.002 on the left-hand side takes 
the time dilation of the i?-meson lifetime t' b = ^tb into account. (In other words, branching 
fractions are Lorentz invariant.) The above formula can be accurately approximated by the 
first term in an expansion in j3 2 , which yields [20] 



l^ ) (E )=r( 1 f\E )-^E 3 



d 1 dT [ u B) 
dE E dE 



+ 0(f3 4 ), 
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as long as E is not too close to the kinematical endpoint (i.e., E < (3_M B /2 2.47 GeV). 
The numerical results for the partial decay rate ^T^iE®) in the rest frame of the T(4S) 
resonance are given in Table 2. 

5.3 Cut on hadronic P + 

Cutting on P + samples the same hadronic phase space as a cut on the charged- lept on energy, 
but with much better efficiency [9, 14]. The phase space P + < Ap with the ideal separator 
Ap = M^/Mb ~ 0.66 GeV contains well over half of all B — > X u l~v events. Here we find 
with the default settings 

p(0) + (pkin(l) + phadr(l)x + ^ m(2) + r ^ adr ( 2 )) 

= [ 53.225 + (4.646 - 11.862) + (0.328 - 0.227)] |K&| 2 ps -1 . 

We see a much better convergence of the power series than in the case of a cut on the charged- 
lepton energy, namely 53.225 — 7.216 — 0.100 when grouping the above numbers according to 
their power counting. Once again, the sum r^ 111 ^ + T^ n ^ = 4.973 is very close to the full 
kinematical correction T^ 111 = 4.959 (in units of |Ktb| 2 ps _1 ). 

Often times it is required to impose an additional cut on the charged-lepton energy, as 
leptons that are too soft are difficult to detect. In Table 3 we list results for both E t > 
and E t > 1.0 GeV. For the ideal cut we find that the prediction is quite precise, as the total 
theoretical uncertainty is only about 6.8%. For comparison, the ideal cut for the lepton energy 
is uncertain by about 50%, but rapidly improving as the energy cut is relaxed. 

5.4 Cut on hadronic invariant mass and q 2 

The most efficient separator for the discrimination of B — > X c l~v events is a cut on the 
invariant mass Mx of the hadronic final state, Mx < Mp> [52, 53]. It has also been argued 
[54] that a cut on q 2 can reduce the shape-function sensitivity, since it avoids the collinear 
region in phase space where P_ ^> P + . In order to optimize signal efficiency and theoretical 
uncertainties, it was suggested in [55] to combine a q 2 cut with a cut on hadronic invariant 
mass. 

The theoretical predictions obtained in [54, 55] were based on a conventional OPE calcu- 
lation, which was assumed to be valid for these cuts. The assessment of the shape-function 
sensitivity was based on convolving the tree-level decay rate with a "tree-level shape function" , 
for which two models (a realistic model similar to the ones considered here, and an unrealistic 
5-function model) were employed. The shape-function sensitivity was then inferred from the 
comparison of the results obtained with the two models. The sensitivity to subleading shape 
functions was not considered, since it was assumed to be very small. Since our formalism 
smoothly interpolates between the "shape-function" and "OPE" regions, and since we include 
radiative corrections as well as power corrections as far as they are known, we can estimate 
the sensitivity of a combined Mx-q 2 cut to the leading and subleading shape functions much 
more accurately. Contrary to [55] , we do not find a significant reduction of the shape-function 
sensitivity when adding the q 2 cut to a cut on hadronic invariant mass. 
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Table 3: Partial decay rate T U (A P , E ) for a cut on the hadronic variable P + < A P and lepton 
energy E\ > E , given in units of |Kife| 2 ps _1 - Predictions are based on the shape-function 
parameter values m b = 4.61 GeV, /4 = 0.2 GeV 2 (top) and m b = 4.55 GeV, /z 2 = 0.3 GeV 2 
(bottom). 
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Table 4: Partial decay rate T U (M , q$, E ) for combined cuts M x < M on hadronic invariant 
mass, q 2 > q$ on leptonic invariant mass, and Ei > E on charged-lepton energy, given in units 
of |T4b| 2 ps _1 - Predictions are based on the shape-function parameter values rrib = 4.61 GeV, 
Hi = 0.2 GeV 2 (top) and m b = 4.55 GeV, /x 2 = 0.3 GeV 2 (bottom). 
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In Table 4 we give results for typical cuts on M x and q 2 , with and without including an 
additional cut on charged-lepton energy. Let us study the contributions for the optimal cut 
M x < M D in detail. We find with the default settings 



= [ 58.541 + (8.027 - 9.048) + (2.100 - 0.318)] |Kb| 2 ps _1 - 



Note the almost perfect (accidental) cancellation of the two terms at order l/m b . The resulting 
power series, 58.541 — 1.022 + 1.782, again exhibits good convergence. As previously, the sum 
pkm(i) _|_ -pkm(2) _ 127 is a good approximation to the exact value T^ m = 9.753 (in units of 
iKifcPps -1 )- The analogous analysis for a combined cut M x < 1.7 GeV and q 2 > 8.0 GeV 2 

r^) _|_ (r^ n(1) + r!r dr(1) ) + (rt in(2) + r^ adr(2) ) 

= [ 25.880 + (4.049 - 6.358) + (1.399 - 0.171)] \V ub \ 2 ps^ , 

which means that the power series is 25.880 — 2.309 + 1.228. Here we have T^ 111 '' 1 -' + T u m ^ = 
5.449, which is close to T^ in = 5.160 (in units of iKfcpps -1 ). 

5.5 Cut on s™ ax and E t 

In [56], the BaBar collaboration employed a cut on both E\ > E and a new kinematical 
variable s^ ax < sq, where the definition for s^ ax involves both hadronic and leptonic variables. 
In the P-meson rest frame, it is 

s max = M 2 + q 2 _ ^ ^ + (g5) 

Rewriting the phase space of this cut in the variables P + , P_, p, we find 

0<P+ < min (M B - 2P , v^) , 
P+ < P_ < min ( — ,M B 



P+<P< min (M B - 2E , P_) , (66) 

where it is understood that if q 2 = (M B — P+)(M B — P ) < (M B — v^o) 2 ' then ^ ne interval 
pmm < < pmax mus t eX cluded from the Pi integration. Here 



1 v + ' ; V 2 2M B J V V 2 2Mfi / 

A summary of our findings is given in Table 5. When compared to the pure charged-lepton 
energy cut in Table 1, the additional cut on sg ax eliminates roughly another 20-30% of events. 
However, the hope is that this cut also reduces the sensitivity to the leading shape function, 
which we expect to be sizable for the pure E\ cut. The uncertainty from subleading shape 
functions, however, is almost unaffected by the s^ ax cut. 
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Table 5: Partial decay rate r u (s , E ) for combined cuts s^ ax < so and E\ > E Q , given in units 
of iKiPps -1 - Predictions are based on the shape-function parameter values m b = 4.61 GeV, 
li\ = 0.2 GeV 2 (top) and m b = 4.55 GeV, /i 2 = 0.3 GeV 2 (bottom). 



so [GeV 2 ] 


E [GeV] 


Mean 


Subl. SF 


Pert. 


Total 


3.5 
3.5 
3.5 
3.5 


1.8 
1.9 
2.0 
2.1 


17.39 
15.86 
13.70 
10.78 


±0.62 
±0.63 
±0.66 
±0.73 


+1.54 
-1.36 
+1.33 
-1.18 
+1.05 
-0.94 
+0.71 
-0.64 


+2.08 
-1.96 
+1.94 
-1.84 
+1.77 
-1.71 
+1.62 
-1.59 


3.5 
3.5 
3.5 
3.5 


1.8 
1.9 
2.0 
2.1 


14.57 
13.18 
11.28 
8.77 


±0.60 
±0.61 
±0.64 
±0.69 


+1.25 
-1.09 
+1.06 
-0.92 
+0.82 
-0.71 
+0.54 
-0.46 


+1.87 
-1.77 
+1.76 
-1.68 
+1.63 
-1.58 
+1.54 
-1.51 



5.6 Eliminating weak annihilation contributions 

In Section 3.5 we have argued that a cut on high q 2 , i.e., q 2 < q 2 , will eliminate the effect of 
weak annihilation and remove the uncertainty associated with this contribution. The cutoff 
q 2 should be small enough to exclude the region around q 2 = m 2 , where this contribution is 
concentrated. It is instructive to assess the "cost" of such an additional cut in terms of the 
loss of efficiency and, more importantly, the behavior of the remaining uncertainties. In order 
to do this, we combine the cut q 2 < (Mb — Ma) 2 with either a cut on P + or on Mx- While 
this particular choice for q 2 still leaves some room to improve the efficiency by increasing q 2 , 
it is not desirable to raise the cut much further, since this would threaten the validity of 
quark-hadron duality. 

The results are summarized in Table 6 and can be compared to the previous "pure" P + 
and Mx cuts in Tables 3 and 4. As an example, let us consider the case P + < 0.65 GeV, 
which is close to the charm threshold. Without the additional q 2 cut we found that the total 
theoretical uncertainty (including the weak annihilation error) is 1^8% . When cutting in 
addition on q 2 < (M B — Mo) 2 , the efficiency decreases by about 20% as expected. However, 
due to the absence of the weak annihilation uncertainty, the overall uncertainty decreases to 
lg 4 %. Therefore both strategies result in comparable relative uncertainties, with a slight favor 
for imposing the additional cut from the theoretical point of view. 

While the small reduction of theoretical errors hardly seems worth the effort of imposing 
the q 2 cut, performing an analysis of the type outlined here and comparing its results with 
those obtained without the additional cut may help to corroborate the expectation that the 
weak annihilation contribution is indeed not much larger than what has been found in [44]. 
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Table 6: Examples of partial decay rates with a cut on q 2 < (Mb — Mo) 2 imposed to eliminate 
the weak annihilation contribution. We consider an additional cut on the hadronic variable 
P+ < (top), or on the hadronic invariant mass M x < M (bottom). As before, decay 
rates are given in units of \V u b\ 2 ps _1 . Predictions are based on the shape-function parameters 
m b = 4.61 GeV and /4 = 0.2 GeV 2 . 



A P [GeV] 


Mean 


Subl. SF 


Pert. 


Total 


0.70 
0.65 
0.60 
0.55 
0.50 


39.96 
37.18 
34.05 
30.61 
26.86 


±1.27 
±1.50 
±1.71 
±1.89 
±1.97 


+2.16 
-2.01 
+1.99 
-1.85 
+1.82 
-1.69 
+1.63 
-1.52 
+1.44 
-1.33 


+2.51 
-2.38 
+2.49 
-2.38 
+2.50 
-2.41 
+2.49 
-2.42 
+2.44 
-2.38 


M [GeV] 


Mean 


Subl. SF 


Pert. 


Total 


M D 
1.70 
1.55 


46.75 
40.70 
33.69 


±0.65 
±1.12 
±1.56 


+2.82 
-2.50 
+2.32 
-2.11 
+1.88 
-1.73 


+2.89 
-2.58 
+2.58 
-2.39 
+2.44 
-2.32 



5.7 Dependence on and shape-function sensitivity 

Non-perturbative hadronic physics enters in our approach via the form of the leading and 
subleading shape functions. The strongest sensitivity by far is to the first moment of the 
leading shape function, which determines the HQET parameter A and with it the 6-quark mass. 
Given that the value of m b = m(,(/x*,/x*) can be determined with good precision from other 
sources (such as moments of the leptonic or hadronic invariant mass spectra in B — > X c l~v 
decays), it is instructive to disentangle this dependence from the sensitivity to higher moments 
or, more generally, to the functional form of the shape functions for fixed m b . 
To explore the dependence on m b we define the exponent 

a (mt ) = ^Li = (^M / (*mt) , m 

d\nm b \ T u J \ m b J 

which means that T u ~ (m b ) a . Table 7 shows the values of this exponent over a wide range of 
values of m b for a variety of experimental cuts. To estimate the sensitivity to the functional 
form we scan over a large set of models for the subleading shape functions, and we also study 
the difference between the results obtained using the exponential or the gaussian ansatz for the 
leading shape function. The corresponding variations are added in quadrature and given as a 
relative change in the corresponding partial decay rates (labeled "Functional Form"). In all 
cases, fi 2 = 0.2 GeV 2 is kept fixed. Because we restrict ourselves to only two functional forms 
for the leading shape function in this study, the resulting sensitivities should be interpreted 
with caution. 
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Table 7: Values of the exponent a(m&) for different kinematical cuts. The parameter (j, 2 = 
0.2 GeV 2 is kept fixed. Also quoted is the sensitivity of the partial decay rates to the functional 
form of the shape functions. See text for explanation. 





m b [GeV] 


4.50 


4.55 


4.60 


4.65 


4.70 


M x < M D 


a 

Functional Form 


9.5 
1.4% 


8.8 
1.1% 


8.2 
0.8% 


7.7 

0.5% 


7.3 
0.4% 


M x < 1.7 GeV 


a 

Functional Form 


12.5 
2.9% 


11.5 
2.6% 


10.5 
2.2% 


9.7 
1.9% 


8.9 
1.6% 


M x < 1.7 GeV 
q 2 > 8 GeV 2 


a 

Functional Form 


10.3 
2.0% 


9.8 
1.7% 


9.3 
1.5% 


9.0 
1.4% 


8.7 
1.4% 


q 2 > {M B - M D ) 2 


a 

Functional Form 


11.4 

5.0% 


11.1 

4.4% 


10.9 
4.0% 


10.8 
3.6% 


10.6 
3.2% 


P+ < M 2 D /M B 


a 

Functional Form 


16.7 

5.3% 


15.0 
4.8% 


13.6 
4.4% 


12.2 
4.0% 


11.1 

3.6% 


Ei > 2.2 GeV 


a 

Functional Form 


22.6 
16.2% 


21.0 
13.1% 


19.7 
11.0% 


18.5 
9.3% 


17.4 
7.9% 



The entries in the table are listed in roughly the order of increasing sensitivity to and 
to the functional form of the shape functions, with the hadronic invariant mass cut showing 
the least sensitivity and the lepton energy cut exhibiting the largest one. To some extent this 
reflects the different efficiencies (or "inclusiveness" ) of the various cuts. It is reassuring that 
a ?s 10 for the pure q 2 cut, in accordance with the findings of [38, 39]. Perhaps somewhat 
surprisingly, for this cut a substantial sensitivity to shape-function effects remains even for 
fixed m b and /i 2 . It is well known that the partial rate with a cut q 2 > (M B — M D ) 2 can 
be calculated using a local OPE in powers of Aq CD /m c [38, 54], thereby avoiding the notion 
of shape-function sensitivity. Differences between the functional forms of the shape functions 
in our approach correspond to effects that are formally of order l/m? c and higher. It is not 
unreasonable that these effects should be of order 3-5%. 

We also checked that for much more relaxed cuts the value of a(rrib) tends to 4.8, as stated 
in (58). For example, for a cut P + < Ap we find (with m& = 4.61 GeV and /i 2 = 0.2 GeV 2 ): 



A P [GeV] 


0.6 0.8 1.0 1.2 1.6 2.0 3.0 M B 


a 


15.4 9.8 7.0 5.8 5.1 5.0 4.9 4.8 



6 Conclusions 

A high-precision measurement of the parameters of the unitarity triangle is an ongoing quest, 
which necessitates the close cooperation of theory and experiment. The determination of 
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\V u b\ from inclusive B — > X u l~v decay requires the measurement of partial decay rates with 
kinematical cuts that eliminate the large background from B — > X c l~v decay, as well as 
theoretical predictions for such quantities. To this end, it is desirable to have a theoretical 
description of the triple differential decay rate, which can be used for predicting arbitrary 
partial rates obtained after integrating over certain regions of phase space. One problem in 
providing such a description is that the power-counting rules of the heavy-quark expansion 
are different in different kinematical domains. In this paper we have overcome this difficulty 

In the shape-function region, our results are in agreement with QCD factorization theorems, 
and perturbative effects have been separated from non-perturbative shape functions. When the 
allowed phase space extends over a large domain, our results smoothly reduce to the expressions 
obtained from the local operator product expansion. We have presented a formalism in which 
event distributions and partial decay rates are expressed without explicit reference to partonic 
quantities such as the 6-quark mass. The sensitivity to such hadronic parameters enters 
indirectly, via the moments of shape functions. The most important non-perturbative object, 
namely the leading-order shape function, can be extracted from the photon spectrum in B — > 
X s 7 decay. This is analogous to extractions of parton distribution functions from fits to data 
on deep inelastic scattering. In this way, the dominant uncertainty from our ignorance about 
bound-state effects in the B meson is turned into an experimental uncertainty, which will 
reduce with increasing accuracy of the experimental data on the photon spectrum. Residual 
hadronic uncertainties are power suppressed in the heavy-quark expansion. 

One goal of this paper was to present a detailed framework in which this program can 
be carried out. We have given formulae that can be readily used for the construction of an 
event generator, as well as to estimate the remaining theoretical uncertainties in a robust and 
automated fashion. 

In practice the leading shape function needs to be parameterized. We have suggested three 
different functional forms, which can be used to fit the data of the B — > X s j photon spectrum. 
Once the data is accurately described by a choice of the shape functions, this function can be 
used in the predictions for partial B — > X u l~v rates and spectra. Subleading shape functions 
give rise to theoretical uncertainties starting at the level of l/m b power corrections. We have 
estimated these uncertainties using a large set of models, each of which obeys the known tree- 
level moment relations, but which are very different in their functional form. A second error 
estimate is determined by the residual renormalization-scale dependence. We also considered 
uncertainties from weak annihilation effects, which in principle can be avoided by cutting away 
the region of phase space in which they contribute. We have suggested a cut on high leptonic 
invariant mass, which accomplishes just that. 

The second half of this paper contains detailed numerical predictions for a variety of partial 
B — > X u l~v decay rates with different kinematical cuts, including cuts on the charged-lepton 
energy (both in the rest frame of the B meson and of the T(4S) resonance), on the hadronic 
quantity P + = Ex — \Px\, on My, on q 2 , and on various combinations of these variables. 
Along with our predictions for the rates we have presented a complete analysis of theoretical 
uncertainties. Once the data on the B — > X s ^ photon spectrum are sufficiently precise to 
accurately determine the leading-order shape function, a determination of \V u b\ with theoretical 
uncertainties at the 5-10% level now seems feasible. 
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A Perturbative Expressions 
A.l Anomalous dimensions 

Here we list the known perturbative expansions of the /3-function and relevant anomalous 
dimensions. We work in the MS scheme and define 



da s (fi) 



din// 

oo 



n=0 



n=0 



Alt, 
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n=0 



47T, 
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as the expansion coefficients for the /^-function, the leading-order SCET current anomalous 
dimension, and the cusp anomalous dimension. To three-loop order, the /5-function reads [57] 
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where n/ = 4 is the number of light flavors, Ca = 3 and GV = 4/3. The three-loop expression 
for the cusp anomalous dimension has recently been obtained in [58]. The coefficients read 
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The SCET anomalous dimension 7 is explicitly known only to one-loop order. However, the 
two-loop coefficient can be extracted by noting that 7 is related to the axial-gauge anomalous 
dimension in deep inelastic scattering [11]. The result is 



7o = -5C F , 
-8CV 
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A. 2 Evolution factor 

The exact expression for the evolution factor reads 

In U (jj, h , Hi) = 2Sr(//h, A*») - 2a r (//ft, hi 



2ay(/// l ,/i i ) , 



(73) 
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where the functions of the right-hand side are solutions to the renormalization-group equations 



-j^— S r (u,/i) = -r cusp (a s (/i)) ln^, 



d 



-a r (v,ii) = -r cusp (a s (/i)) , 



-a y (i/,/x) = -V(a s (/i)) , 



am/x am/x 

with boundary conditions Siy, /x) = etc. at /x = v. These equations can be integrated using 
that d/aTn/x = (3(a s ) d/ da s . The solutions are 



SrM = -J da E^MJ 
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and similarly for ay. 

Next, we give explicit results for the Sudakov exponent Sr and the functions ar and a 7 
in (73) at next-to-leading order in renormalization-group improved perturbation theory. We 
obtain 
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To Po, 

and similarly for a 7 . The next-to-leading order expressions for the Sudakov exponent Sr 
contains the three-loop coefficients (3 2 and T 2 . With r = a s (ji) / a s (u), it reads 
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The next-to-leading-logarithmic evolution factor £/(/i/j,/Xj) can be obtained by combining the 
above expressions according to (73) and expanding out terms of order a s . 



B Partially integrated decay rates 

With the exception of the combined cut on the lepton energy Ei and the hadronic quantity 
s^ ax studied in Section 5.5, all other partial rates investigated in our analysis can be derived 
by first integrating the triple differential decay rate (23) over the lepton energy Ei > E and 
P_ < p max analytically, where the quantity P max (and in principle even E ) may depend on 
the value of P + . The remaining integration over P + is then performed numerically. In such a 
situation, we need to evaluate the partially integrated decay rate 

li = L dP -L dP 'mi^- (78) 



41 



Changing variables from P_ to y defined in (25), the constraint P_ < P™ ax translates into the 
integration domain < y < y max , where in analogy to (25) we define 

_ - p+ _ pr x - p+ _ 1 2p 

y™*~ M B -P + ' M B -P + " m b -p + - { > 

From the phase-space relation (2) it follows that a cut on the lepton energy has no effect if 
yo > 2/max- The result of performing the integrations in (78) can be written as 

max 

^ ' " \r^(y ) + rf ; y max >y , 

where 

rftl/i) = (M B - P + ) 5 Etyx fc , ^) £dyy 2 - 2ar [(3 - 2y) ^ + 6(1 - y) % + y <F 3 ] , 

= S^!! {Mb - P + ) 5 ^) f^dyy-^yo 

x [(63/(1 + yo) - Qy 2 - J/o(3 + 2y )) + 6y(l - y) 5F 2 + y (3y - 2y ) J 3 ] • (81) 

When the kinematical power corrections in (30) are expanded as in (32) and (33), the resulting 
integrals over y can be expressed in terms of the master functions I n (b, z) given in eq. (86) of 
[9]. The resulting expressions are used to obtain the numbers in the various tables in Section 5. 

We now list the values of yo an d y m ax for the different cuts studied in Section 5. Whenever 
a cut Ei > Eq on the charged-lepton energy is applied, we have 

^~jBp + - (82) 

For an additional cut P + < A P , we have y max = 1 and < P + < min(A P , M B — 2E ). For a 
cut on hadronic invariant mass, M x < M , we have 

ymax M B -P+ K ' 

and < P + < min(M , M B — 2E ). For a cut on leptonic invariant mass, q 2 > q%, we have 

Vmax = 1 " (M B % P + ) 2 (84) 

and < P + < min(M B — q , M B — 2E ). Finally, for the combined M x -q 2 cut we take the 
minimum of the previous two y max values. 
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